PARETO OPTIMA OF 
MULTICRITERIA INTEGER LINEAR PROGRAMS 



JESUS A. DE LOERA, RAYMOND HEMMECKE, AND MATTHIAS KOPPE 

Abstract. We settle the computational complexity of fundamental 
questions related to multicriteria integer linear programs, when the di- 
mensions of the strategy space and of the outcome space are considered 
fixed constants. In particular we construct: 

1. polynomial-time algorithms to exactly determine the number of 
Pareto optima and Pareto strategies; 

2. a polynomial-space polynomial-delay prescribed-order enumera- 
tion algorithm for arbitrary projections of the Pareto set; 

3. an algorithm to minimize the distance of a Pareto optimum from a 
prescribed comparison point with respect to arbitrary polyhedral 
norms; 

4. a fully polynomial-time approximation scheme for the problem of 
minimizing the distance of a Pareto optimum from a prescribed 
comparison point with respect to the Euclidean norm. 



1. Introduction 

Let A = (ciij) be an integral m x re-matrix and b £ Z m such that the 
convex polyhedron P = { u € W 1 : Au < b } is bounded. Given k linear 
functional /i,/2,...,/t £ Z", we consider the multicriterion integer linear 
programming problem 

vmin (/i(u),/ 2 (u),...,/fc(u)) 
subject to Au < b (1) 

uez n 

where vmin is defined as the problem of finding all Pareto optima and a 
corresponding Pareto strategy. For a lattice point u the vector f(u) = 
(/i(u), . . . , /fe(u)) is called an outcome vector. Such an outcome vector is a 
Pareto optimum for the above problem if and only if there is no other point u 
in the feasible set such that /i(u) < fi(u) for all i and fj(u) < /j(u) for 
at least one index j. The corresponding feasible point u is called a Pareto 
strategy. Thus a feasible vector is a Pareto strategy if no feasible vector 
can decrease some criterion without causing a simultaneous increase in at 
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least one other criterion. For general information about the multicriteria 
problems see, e.g., [8, 12]. 

In general multiobjective problems the number of Pareto optimal solu- 
tions may be infinite, but in our situation the number of Pareto optima 
and strategies is finite. There are several well-known techniques to gen- 
erate Pareto optima. Some popular methods used to solve such problems 
include, e.g., weighting the objectives or using a so-called global criterion 
approach (see [6]). In abnormally nice situations, such as multicriteria lin- 
ear programs [9], one knows a way to generate all Pareto optima, but most 
techniques reach only some of the Pareto optima. 

The purpose of this article is to study the sets of all Pareto optima and 
strategies of a multicriterion integer linear program using the algebraic struc- 
tures of generating functions. The set of Pareto points can be described as 
the formal sum of monomials 

J2{ z v : u G P n Z n and v = f(u) G Z k is a Pareto optimum }. (2) 

Our main theoretical result states that, under the assumption that the num- 
ber of variables is fixed, we can compute in polynomial time a compact ex- 
pression for the huge polynomial above, thus all its Pareto optima can in 
fact be counted exactly. The same can be done for the corresponding Pareto 
strategies when written in the form 

^{ x u : u£PnZ" and f (u) is a Pareto optimum } . (3) 

Theorem 1. Let A G Z mxn ; a d-vector h, and linear functions fi, . . . , ft G 
Z n be given. There are algorithms to perform the following tasks: 

(i) Compute the generating function (2) of all the Pareto optima as a sum 
of rational functions. In particular we can count how many Pareto 
optima are there. If we assume k and n are fixed, the algorithm runs 
in time polynomial in the size of the input data. 

(ii) Compute the generating function (3) of all the Pareto strategies as a 
sum of rational functions. In particular we can count how many Pareto 
strategies are there in P. If we assume k and n are fixed, the algorithm 
runs in time polynomial in the size of the input data. 

(iii) Generate the full sequence of Pareto optima ordered lexicographically 
or by any other term ordering. If we assume k and n are fixed, the 
algorithm runs in polynomial time on the input size and the num- 
ber of Pareto optima. (More strongly, there exists a polynomial-space 
polynomial-delay prescribed-order enumeration algorithm. ) 

In contrast it is known that for non-fixed dimension it is #P-hard to 
enumerate Pareto optima and NP-hard to find them [7, 13]. The proof of 
Theorem 1 parts (i) and (ii) will be given in section 2. It is based on the 
theory of rational generating functions as developed in [1, 2]. Part (iii) of 
Theorem 1 will be proved in section 3. 



MULTICRITERIA INTEGER LINEAR PROGRAMS 



3 



For a user that knows some or all of the Pareto optima or strategies, a goal 
is to select the "best" member of the family. One is interested in selecting one 
Pareto optimum that realizes the "best" compromise between the individual 
objective functions. The quality of the compromise is often measured by the 
distance of a Pareto optimum v from a user-defined comparison point v. For 
example, often users take as a good comparison point the so-called ideal point 
v ideai g jk Q f multicriterion problem, which is defined as 

^ dcal = min{ /i(u) : u G P n 1 n }. 

The criteria of comparison with the point v are quite diverse, but some 
popular ones include computing the minimum over the possible sums of 
absolute differences of the individual objective functions, evaluated at the 
different Pareto strategies, from the comparison point v, i.e., 

/(u) = |A(u) - + • • • + |/ fc (u) - v k \, (4a) 

or the maximum of the absolute differences, 

f(u) =max{|/i(u) -t>i|,...,|/fc(u) -v k \}, (4b) 

over all Pareto optima (/i(u), . . . , f k (u)). Another popular criterion, some- 
times called the global criterion, is to minimize the sum of relative distances 
of the individual objectives from their known minimal values, i.e., 

/<») = *^ + -" + *^- (-) 

\ v i I \ v k I 

We stress that if we take any one of these functions as an objective func- 
tion of an integer program, the optimal solution will be a non-Pareto solution 
of the multicriterion problem (1) in general. In contrast, we show in this 
paper that by encoding Pareto optima and strategies as a rational function 
we avoid this problem, since we evaluate the objective functions directly on 
the space of Pareto optima. 

All of the above criteria (4) measure the distance from a prescribed point 
with respect to a polyhedral norm. In section 4, we prove: 

Theorem 2. Let the dimension n and the number k of objective functions 
be fixed. Let a multicriterion integer linear program (1) be given. Let a 
polyhedral norm \\-\\q be given by the vertex or inequality description of its 
unit ball QCR fc . Finally, let a prescribed point v € T, k be given. 

(i) There exists a polynomial-time algorithm to find a Pareto optimum v 
of (1) that minimizes the distance ||v — v||q from the prescribed point. 

(ii) There exists a polynomial- space polynomial- delay enumeration algo- 
rithm for enumerating the Pareto optima of (1) in the order of in- 
creasing distances from the prescribed point v. 

Often users are actually interested in finding a Pareto optimum that min- 
imizes the Euclidean distance from a prescribed comparison point v, 

/(u) = ^\f 1 (u)-v l \ 2 + --- + \f k (u)-v k \ 2 , (5) 
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but to our knowledge no method of the literature gives a satisfactory solution 
to that problem. In section 4, however, we prove the following theorem, 
which gives a very strong approximation result. 

Theorem 3. Let the dimension n and the number k of objective functions 
be fixed. There exists a fully polynomial-time approximation scheme for the 
problem of minimizing the Euclidean distance of a Pareto optimum of (1) 
from a prescribed comparison point v£Z fc . 

We actually prove this theorem in a somewhat more general setting, us- 
ing an arbitrary norm whose unit ball is representable by a homogeneous 
polynomial inequality. 

2. The rational function encoding of all Pareto optima 

We give a very brief overview of the theory of rational generating functions 
necessary to establish Theorem 1. For full details we recommend [1, 2, 3, 5] 
and the references therein. In 1994 Barvinok gave an algorithm for counting 
the lattice points in P = { u 6 1R™ : Au < b } in polynomial time when 
the dimension n of the feasible polyhedron is a constant [1] . The input for 
Barvinok's algorithm is the binary encoding of the integers aij and bi, and 
the output is a formula for the multivariate generating function 

g(P;x) = Y, x " 

where x u is an abbreviation of x^x^ 2 ■ ■ -x^™- This long polynomial with 
exponentially many monomials is encoded as a much shorter sum of rational 
functions of the form 

x c * 

^ P > X ) = 7i (l- X d*i)(l_xd ! 2)...(l_ x d ln )- (6) 

Barvinok and Woods in 2003 further developed a set of powerful manipula- 
tion rules for using these short rational functions in Boolean constructions 
on various sets of lattice points. 

Throughout the paper we assume that the polyhedron P = { u € M. n : 
Au < b } is bounded. We begin by recalling some useful results of Barvinok 
and Woods (2003): 

Theorem 4 (Intersection Lemma; Theorem 3.6 in [3]). Let £ be a fixed in- 
teger. Let Si, S2 be finite subsets ofZ n . Let <?(Si;x) and g(S2] x) be their 
generating functions, given as short rational functions with at most I bino- 
mials in each denominator. Then there exists a polynomial time algorithm, 
which computes 



g(Si n5 2 ;x) = ^ 



li- 



(1 -x d *i) ... (1 - X d ") 

with s < 2£, where the ji are rational numbers, Cj,dy are nonzero integer 
vectors, and I is a polynomial- size index set. 
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The following theorem was proved by Barvinok and Woods using Theo- 
rem 4: 

Theorem 5 (Boolean Operations Lemma; Corollary 3.7 in [3]). Let m and 

t be fixed integers. Let S±, S2, • • • , S m be finite subsets ofZ n . Let g(Sf, x) for 
i = 1, . . . , m be their generating functions, given as short rational functions 
with at most I binomials in each denominator. Let a set S C Z n be defined 
as a Boolean combination of S\, . . . , S m (i.e., using any of the operations U, 
n, \). Then there exists a polynomial time algorithm, which computes 



where s = s(£,m) is a constant, the ji are rational numbers, Ci,djj are 
nonzero integer vectors, and L is a polynomial-size index set. 

We will use the Intersection Lemma and the Boolean Operations Lemma 
to extract special monomials present in the expansion of a generating func- 
tion. The essential step in the intersection algorithm is the use of the 
Hadamard product [3, Definition 3.2] and a special monomial substitution. 
The Hadamard product is a bilinear operation on rational functions (we 
denote it by *). The computation is carried out for pairs of summands as 
in (6). Note that the Hadamard product rrt\ * 777-2 of two monomials mi, 7772 
is zero unless mi = mi- 

Another key subroutine introduced by Barvinok and Woods is the follow- 
ing Projection Theorem. 

Theorem 6 (Projection Theorem; Theorem 1.7 in [3]). Assume the dimen- 
sion n is a fixed constant. Consider a rational polytope P C W 1 and a linear 
map T: Z n -» Z k . There is a polynomial time algorithm which computes a 
short representation of the generating function f(T(P n Z n );x). 

One has to be careful when using earlier Lemmas (especially the Projec- 
tion Theorem) that the sets in question are finite. The proof of Theorem 1 
will require us to project and intersect sets of lattice points represented by 
rational functions. We cannot, in principle, do those operations for infinite 
sets of lattice points. Fortunately, in our setting it is possible to restrict our 
attention to finite sets. 

Finally, one important comment. If we want to count the points of a 
lattice point set S, such as the set of Pareto optima, it would apparently 
suffice to substitute 1 for all the variables xi of the generating function 



to get the specialization \S\ = g(S; x = 1). But this cannot be done directly 
due to the singularities in the rational function representation. Instead, 





6 



J. A. DE LOERA, R. HEMMECKE, AND M. KOPPE 



choose a generic vector A = (Ai, . . . , X n ) and substitute each of the vari- 
ables Xi by e tXi . Then we get 

9(S,e ) = Yl^Jl _ e t(A,d il >)(i _ e *<^.di2>) ... (1 - e *( A ' d -)) ' 

Counting the number of lattice points is the same as computing the constant 
terms of the Laurent series for each summand and adding them up. This 
can be done using elementary complex residue techniques (see [2]). 

Proof of Theorem 1, part (i) and (ii). The proof of part (i) has three steps: 
Step 1. For i = 1, . . . , k let V{ € Z be an upper bound of polynomial 
encoding size for the value of fi over P. Such a bound exists because of 
the boundedness of P, and it can be computed in polynomial time by linear 
programming. We will denote the vector of upper bounds by v £ We 
consider the truncated multi- epigraph of the objective functions /i, •••,/& 
over the linear relaxation of the feasible region P, 

pf . = { ( u , v) € R n x R k : u e P, 

Vi > Vi > /i(u) for i = 1, . . . , k }, 

which is a rational convex polytope in R n x R k . Let V- C Z k denote the 
integer projection of Pjr j, on the v variables, i.e., the set 

V- = {v £Z k :3ueZ n with (u, v) € Pfc f fl {Z n x Z fc ) } . (8) 

Clearly, the vectors in V- are all integer vectors in the outcome space which 
are weakly dominated by some outcome vector (/i(u), /2(u), . . . , /fc(u)) for 
a feasible solution x in PnZ n ; however, we have truncated away all outcome 
vectors which weakly dominate the computed bound v. Let us consider the 
generating function of V-, the multivariate polynomial 

g(V^;z) =^{z v :vG^}. 

In the terminology of polynomial ideals, the monomials in g(V-;z) form a 
truncated ideal generated by the Pareto optima. By the Projection Theorem 
(our Theorem 6), we can compute g(V-;z) in the form of a polynomial-size 
rational function in polynomial time. 

Step 2. Let y Parcto c Z fc denote the set of Pareto optima. Clearly we have 

y Parcto = (y> \ ^ + y>^ f| • • • Pi (V~ \ (efc + V~)) , 

where € 1> k denotes the i-th unit vector and 

e i + V^ = {e i + v:veV r ^}. 
The generating function g,(y Pareto ; z ) can De computed by the Boolean Op- 
erations Lemma (Theorem 5) in polynomial time from g(V-;z) as 

5(y Parcto. z) = ( 9( y> ;z) - 5 (y>; Z ) * Zl g(V^,z)) 

* — *(g{V*;z)-g(yZ;z)*z k g(y*;z)) 1 
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where * denotes taking the Hadamard product of the rational functions. 

Step 3. To obtain the number of Pareto optima, we compute the special- 
ization g(V Parcto ;z = l). This is possible in polynomial time using residue 
techniques as outlined before the beginning of the proof. 

Proof of part (ii). Now we recover the Pareto strategies that gave rise to the 
Pareto optima, i.e., we compute a generating function for the set 

^Pareto = | u G p : u e p n Z" and f(u) is a Pareto optimum }. 

To this end, we first compute the generating function for the set 

^Pareto = { ( Uj v ) e z n x Z fc : v is a p are t point with Pareto strategy u } . 

For this purpose, we consider the multi-graph of the objective functions 
fi,...,fk over P, 

P= = { ( u , v) € R n x R k : u G P, 

lV ' ' (10) 
Vi = fi(u) for i = 1, . . . ,k}. 

Using Barvinok's theorem, we can compute in polynomial time the generat- 
ing function for the integer points in P, 

9 (P;x) = ^{x u :u(;PnP}, 

and also, using the monomial substitution Xj — > xjzf 1 ^^ ■ ■ ■ zl k ^ ej ^ for all j, 
the generating function is transformed into 

fr*-*> = E( xUzV : ( u ' v ) e P h,..J k n ^ l >< ^) h 

where the variables x carry on the monomial exponents the information of 
the u-coordinates of Pj x j k and the z variables of the generating function 
carry the v-coordinates of lattice points in Pjr . Now 

5 (S Parcto ;x,z) = ( 5 (P;x) 5 (y Parcto ;z)) * ^ ... Jfc ; x, z), (11) 

which can be computed in polynomial time for fixed dimension by the the- 
orems outlined early on this section. Finally, to obtain the generating func- 
tion g([/ Pareto ; x) of the Pareto strategies, we need to compute the projection 
of 5 Pareto into the space of the strategy variables u. Since the projection is 
one-to-one, it suffices to compute the specialization 

5 (C/ Parcto ;x)= 5 (5 Parcto ;x,z = l), 

which can be done in polynomial time. □ 

3. Efficiently listing all Pareto optima 

The Pareto optimum that corresponds to the "best" compromise between 
the individual objective functions is often chosen in an interactive mode, 
where a visualization of the Pareto optima is presented to the user, who 
then chooses a Pareto optimum. Since the outcome space frequently is of a 
too large dimension for visualization, an important task is to list (explicitly 
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enumerate) the elements of the projection of the Pareto set into some lower- 
dimensional linear space. 

It is clear that the set of Pareto optima (and thus also any projection) is 
of exponential size in general, ruling out the existence of a polynomial-time 
enumeration algorithm. In order to analyze the running time of an enumer- 
ation algorithm, we must turn to output-sensitive complexity analysis. 

Various notions of output-sensitive efficiency have appeared in the lit- 
erature; we follow the discussion of [10]. Let W C IP be a finite set to 
be enumerated. An enumeration algorithm is said to run in polynomial 
total time if its running time is bounded by a polynomial in the encod- 
ing size of the input and the output. A stronger notion is that of in- 
cremental polynomial time: Such an algorithm receives a list of solutions 
Wi,...,wjv £ W as an additional input. In polynomial time, it outputs 
one solution w € W \ {wi, . . . , wyv} or asserts that there are no more so- 
lutions. An even stronger notion is that of a polynomial- delay algorithm, 
which takes only polynomial time (in the encoding size of the input) before 
the first solution is output, between successive outputs of solutions, and af- 
ter the last solution is output to the termination of the algorithm. Since 
the algorithm could take exponential time to output all solutions, it could 
also build exponential-size data structures in the course of the enumera- 
tion. This observation gives rise to an even stronger notion of efficiency, a 
polynomial- space polynomial- delay enumeration algorithm. 

We also wish to prescribe an order, like the lexicographic order, in which 
the elements are to be enumerated. We consider term orders -<r on mono- 
mials y w that are defined as in [11] by a non- negative integral p x p-matrix 
R of full rank. Two monomials satisfy y Wl -<_r y W2 if and only if Rw\ 
is lexicographically smaller than Rw2- I n other words, if ri, . . . ,r n denote 
the rows of R, there is some j € {1, ...,n} such that (rj,wi) = (rj,W2) 
for i < j, and (r^Wi) < (rj,w 2 ). For example, the unit matrix R = I n 
describes the lexicographic term ordering. 

We prove the existence of a polynomial-space polynomial-delay prescribed- 
order enumeration algorithm in a general setting, where the set W to be enu- 
merated is given as the projection of a set presented by a rational generating 
function. 

Theorem 7. Let the dimension k and the maximum number i of binomials 
in the denominator be fixed. 

Let V C Z fc be a bounded set of lattice points with V C [— M, M] k , given 
only by the bound M € Z + and its rational generating function encoding 
g(V; z) with at most I binomials in each denominator. Let 

W = { w € IT : 3t € 1 k ~ p such that (t, w) £ V } 

denote the projection of V onto the last p components. Let -<r be the term 
order on monomials in y\, . . . , y p induced by a given matrix R G N pxp . 

There exists a polynomial- space polynomial- delay enumeration algorithm 
for the points in the projection W , which outputs the points of W in the 
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order given by The algorithm can be implemented without using the 

Projection Lemma. 

We remark that Theorem 7 is a stronger result than what can be ob- 
tained by the repeated application of the monomial-extraction technique of 
Lemma 7 from [4], which would only give an incremental polynomial time 
enumeration algorithm. 

Proof. We give a simple recursive algorithm that is based on the iterative 
bisection of intervals. 

Input: Lower and upper bound vectors 1, u € Z p . 

Output: All vectors winlf with 1 < Rw < u, sorted in the 

order 

1. If the set W n { w : 1 < Rw < u } is empty, do nothing. 

2. Otherwise, if 1 = u, compute the unique point w € Z fc 
with Rw = 1 = u and output w. 

3. Otherwise, let j be the smallest index with lj ^ Uj. 
We bisect the integer interval {lj, . . . , Uj} evenly into 

{lj, . . . , nij} and {rrij + 1, . . . , Uj}, where rrij = . 

We invoke the algorithm recursively on the first part, 
then on the second part, using the corresponding lower 
and upper bound vectors. 

We first need to compute appropriate lower and upper bound vectors 1, u 
to start the algorithm. To this end, let iV be the largest number in the 
matrix R and let 1 = -pMNl and u = pMNl. Then 1 < Rw < u holds for 
all w £ W. Clearly the encoding length of 1 and u is bounded polynomially 
in the input data. 

In step 1 of the algorithm, to determine whether 

l^n{w:l<i?w<u} = 0, (12) 

we consider the polytope 

Q l u = [-M, M} k ~ p x { w € M p : 1 < Rw < u } C E fc , (13) 

a parallelelepiped in M. k . Since W is the projection of V and since V C 
[-M,M] k , we have (12) if and only if V Cl Qi, u = 0. The rational gen- 
erating function g(Qi tU ] z ) can be computed in polynomial time. By us- 
ing the Intersection Lemma, we can compute the rational generating func- 
tion g(V fl Qi )U ; z) in polynomial time. The specialization g(V n Qi, u ; z = 1) 
can also be computed in polynomial time. It gives the number of lattice 
points in V n Qi, u ; in particular, we can decide whether V n Q\ tU = 0. 

It is clear that the algorithm outputs the elements of W in the order given 
by -<r- We next show that the algorithm is a polynomial-space polynomial- 
delay enumeration algorithm. The subproblem in step 1 only depends on 
the input data as stated in the theorem and on the vectors 1 and u, whose 
encoding length only decreases in recursive invocations. Therefore each of 
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the subproblems can be solved in polynomial time (thus also in polynomial 
space) . 

The recursion of the algorithm corresponds to a binary tree whose nodes 
are labeled by the bound vectors 1 and u. There are two types of leaves in 
the tree, one corresponding to the "empty-box" situation (12) in step 1, and 
one corresponding to the "solution-output" situation in step 2. Inner nodes 
of the tree correspond to the recursive invocation of the algorithm in step 3. 
It is clear that the depth of the recursion is 0(plog(pMN)), because the 
integer intervals are bisected evenly. Thus the stack space of the algorithm 
is polynomially bounded. Since the algorithm does not maintain any global 
data structures, the whole algorithm uses polynomial space only. 

Let Wj € W be an arbitrary solution and let Wj+i be its direct successor 
in the order We shall show that the algorithm only spends polynomial 
time between the output of Wj and the output of w, + i. The key property 
of the recursion tree of the algorithm is the following: 

Every inner node is the root of a subtree that contains at , 
least one solution-output leaf. 

The reason for that property is the test for situation (12) in step 1 of the 
algorithm. Therefore, the algorithm can visit only 0(p\og(pMN)) inner 
nodes and empty-box leaves between the solution-output leaves for w; and 
Wj-t-i. For the same reason, also the time before the first solution is output 
and the time after the last solution is output are polynomially bounded. □ 

The following corollary, which is a stronger formulation of Theorem 1 (iii) , 
is immediate. 

Corollary 8. Let n and k be fixed integers. There exist polynomial- space 
polynomial- delay enumeration algorithms to enumerate the set of Pareto op- 
tima of the multicriterion integer linear program (1), the set of Pareto strate- 
gies, or arbitrary projections thereof in lexicographic order ( or an arbitrary 
term order). 

Remark 9. We remark that Theorem 7 is of general interest. For in- 
stance, it also implies the existence of a polynomial-space polynomial-delay 
prescribed-order enumeration algorithm for Hilbert bases of rational poly- 
hedral cones in fixed dimension. 

Indeed, fix the dimension d and let C = cone{bi, . . . , b n } C M. d be a 
pointed rational polyhedral cone. The Hilbert basis of C is defined as the 
inclusion-minimal set H C C n Z d which generates C H 7L d as a monoid. For 
simplicial cones C (where bi,...,b n are linearly independent), Barvinok 
and Woods [3] proved that one can compute the rational generating function 
g(H; z) (having a constant number of binomials in the denominators) of the 
Hilbert basis of C n Z d using the Projection Theorem. The same technique 
works for non-simplicial pointed cones. Now Theorem 7 gives a polynomial- 
space polynomial-delay prescribed-order enumeration algorithm. 



MULTICRITERIA INTEGER LINEAR PROGRAMS 



11 



4. Selecting a Pareto optimum using global criteria 

Now that we know that all Pareto optima of a multicriteria integer linear 
programs can be encoded in a rational generating function, and that they can 
be listed efficiently on the output size, we can aim to apply selection criteria 
stated by a user. The advantage of our setup is that when we optimize a 
global objective function it guarantees to return a Pareto optimum, because 
we evaluate the global criterion only on the Pareto optima. Let us start 
with the simplest global criterion which generalizes the use of the l\ norm 
distance function: 

Theorem 10. Let the dimension k and the maximum number £ of binomials 
in the denominator be fixed. 

Let V C Z fc be a bounded set of lattice points with V C [— M, M] n+k , given 
only by the bound M G Z + and its rational generating function encoding 
g(V; z) with at most t binomials in the denominators. 

Let Q C W k be a rational convex central- symmetric polytope with G 
int Q, given by its vertex or inequality description. Let the polyhedral norm 
\\-\\q be defined using the Minkowski functional 



Finally, let a prescribed point v G Z fc be given. 

(i) There exists a polynomial-time algorithm to find a point v G V that 
minimizes the distance dg(v, v) = ||v — v||q from the prescribed point. 

(ii) There exists a polynomial- space polynomial- delay enumeration algo- 
rithm for enumerating the points of V in the order of increasing dis- 
tances dQ from the prescribed point v, refined by an arbitrary term 
order given by a matrix R G N kxk . 

Theorem 2, as stated in the introduction, is an immediate corollary of 
this theorem. 

Proof. Since the dimension k is fixed, we can compute an inequality descrip- 
tion 



of Q with A £ ^ mxfc anc i b G Z fc in polynomial time, if Q is not already 
given by an inequality description. Let v € V be arbitrary; then 



y|| o =inf{A>0:y G XQ}. 



(15) 



Q = {y GR k :Ay <b} 



dg(v,v) = ||v- v|| Q 

= inf { A>0:v-vGAQ} 

= min{ A > : Ab > A(v - v) }. 

Thus there exists an index i G {1, . . . , m} such that 
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so (2q(v,v) is an integer multiple of 1/bi. Hence for every v 6 V, we have 
that 

d Q(v,v)e-— -i TT Z +' ( 16 ) 

lcm(6i,...,6 m ) 

where lcm(&i, . . . ,b m ) clearly is a number of polynomial encoding size. On 
the other hand, every v £ V certainly satisfies 

dg(v, v) < ka(M + max{|£>i| , . . . , \vd\}) (17) 

where a is the largest number in A, which is also a bound of polynomial 
encoding size. 

Using Barvinok's algorithm, we can compute the rational generating func- 
tion g(v + XQ; z) for any rational A of polynomial enoding size in polyno- 
mial time. We can also compute the rational generating function g(V D 
(v + XQ);z) using the Intersection Lemma. By computing the specializa- 
tion g(V n (v + XQ);z = 1), we can compute the number of points in 

V PI (v + XQ), thus we can decide whether this set is empty or not. 
Hence we can employ binary search for the smallest A > such that 

V n (v + XQ) is nonempty. Because of (16) and (17), it runs in polynomial 
time. By using the recursive bisection algorithm of Theorem 7, it is then 
possible to construct one Pareto optimum in V n (v + XQ) for part (i), or to 
construct a sequence of Pareto optima in the desired order for part (ii). □ 

Now we consider a global criterion using a distance function corresponding 
to a non-polyhedral norm like the Euclidean norm ||-|| 2 (or any other £ p -norm 
for 1 < p < oo). We are able to prove a very strong type of approximation 
result, a so-called fully polynomial-time approximation scheme (FPTAS), in 
a somewhat more general setting. 

Definition 11 (FPTAS). Consider the optimization problems 

max{/(v) : v G V}, (18a) 
min{/(v) : v € V}. (18b) 

A fully polynomial-time approximation scheme (FPTAS) for the maximiza- 
tion problem (18a) or the minimization problem (18b), respectively, is a 
family {^4 e :e€Q,e>0}of approximation algorithms A e , each of which 
returns an e- approximation, i.e., a solution v e € V with 

/(v e ) > (1 - e)/* where /* = max/(v), (19a) 

vgV 

or, respectively, 

/(v e ) < (1 + e)/* where /*=min/(v), (19b) 

vev 

such that the algorithms A e run in time polynomial in the input size and \. 
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Figure 1 . A set defining a pseudo-norm with the inscribed 
and circumscribed cubes aB^ and /3-E>oo (dashed). 

Remark 12. An FPTAS is based on the notion of e-approximation (19), 
which gives an approximation guarantee relative to the value /* of an op- 
timal solution. It is clear that this notion is most useful for objective func- 
tions / that are non-negative on the feasible region V. Since the approxi- 
mation quality of a solution changes when the objective function is changed 
by an additive constant, it is non-trivial to convert an FPTAS for a maxi- 
mization problem to an FPTAS for a minimization problem. 

We shall present an FPTAS for the problem of minimizing the distance of 
a Pareto optimum from a prescribed outcome vector v G Z fc . We consider 
distances d(v, •) induced by a pseudo-norm ||-||g via 

d(v,v) = ||v- v|| Q (20a) 

To this end, let Q C M fc be a compact basic semialgebraic set with G int Q, 
which is described by one polynomial inequality, 

Q={y£R k :q(y)<l}, (20b) 

where q G Q[yi, ■ ■ ■ ,Uk] is a homogeneous polynomial of (even) degree D. 
The pseudo-norm ||-||q is now defined using the Minkowski functional 

||y||g = inf{ A > : y G AQ } (20c) 

Note that we do not make any assumptions of convexity of Q, which would 
make ||-||q a norm. Since Q is compact and G intQ, there exist positive 
rational numbers (norm equivalence constants) a, (3 with 

aB^ CQC fiBoo where = { y G R fc : fly^ < 1 }; (21) 

see Figure 1. 

Now we can formulate our main theorem, which has Theorem 3, which 
we stated in the introduction, as an immediate corollary. 

Theorem 13. Let the dimension n and the number k of objective functions 
be fixed. Moreover, let a degree D and two rational numbers < a < (3 
be fixed. Then there exists a fully polynomial-time approximation scheme 
for the problem of minimizing the distance dq(v, v), defined via (20) by a 
homogeneous polynomial q G Q[yi, • • • , yk] of degree D satisfying (21), whose 
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coefficients are encoded in binary and whose exponent vectors are encoded in 
unary, of a Pareto optimum of (1) from a prescribed outcome vector v G Z fc . 

The proof is based on the following result, which is a more general for- 
mulation of Theorem 1.1 from [5]. 

Theorem 14 (FPTAS for maximizing non-negative polynomials over finite 
lattice point sets). For all fixed integers k (dimension) and s (maximum 
number of binomials in the denominator), there exists an algorithm with 
running time polynomial in the encoding size of the problem and - for the 
following problem. 

Input: Let V C Z fc be a finite set, given by a rational generating function 
in the form 

. a (l-x d *i)...(l-x d -i) 

where the the numbers Sj of binomials in the denominators are at most s. 
Furthermore, let two vectors vl, vjj G Z fc be given such that V is contained 
in the box { v : vl < v < vtj } . 

Let f G Q[i>i, • • • , Ufc] be a polynomial with rational coefficients that is 
non-negative on V , given by a list of its monomials, whose coefficients are 
encoded in binary and whose exponents are encoded in unary. 

Finally, let e G Q. 

Output: Compute a point v e G V that satisfies 

/(v e ) > (1 - e)f* where f* = max/(v). 

In [5] the above result was stated and proved only for sets V that consist 
of the lattice points of a rational polytope; however, the same proof yields 
the result above. 



Proof of Theorem 13. Using Theorem 1, we first compute the rational gener- 
ating function g(y Pareto ; z ) Q f the Pareto optima. With binary search using 
the Intersection Lemma with generating functions of cubes as in section 3, 
we can find the smallest non-negative integer 7 such that 

(v + 7 5 oo )ny Paict V0. (22) 

If 7 = 0, then the prescribed outcome vector v itself is a Pareto optimum, 
so it is the optimal solution to the problem. 

Otherwise, let vo be an arbitrary outcome vector in (v + 7-B00) H y Pareto . 
Then 

7 > ||v - vH^ = inf{ A : v - v G XB^ } 

> inf{ A : v - v G X^Q } = a ||v - v\\ Q , 
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thus ||vo — v||g < y/a. Let 5 = (3y/a. Then, for every vi G M k with 
|| Vi — v || ^ > 5 we have 

5 < || vi - vjl^ = inf{ A : vi - v G XBoo } 

< inf { A : vi - v G A^Q } = /3 ||vi - v|| Q , 

thus 

||vi - v|| Q > 5/(3 = y/a > ||v - v\\ Q ■ 

Therefore, a Pareto optimum v* G y Pareto minimizing the distance d,Q from 
the prescribed outcome vector v is contained in the cube v+5-Boo. Moreover, 
for all points v G v + SBqq we have 

II v - v||q < S/a = (3y/a 2 . 

We define a function / by 

/(v) = (/? 7 /a 2 ) Z) -||v-v||g, (23) 

which is non-negative over the cube v + SB^. Since q is a homogeneous 
polynomial of degree D, we obtain 

/(v) = (/? 7 /a 2 ) D -g(v-v) (24) 

so / is a polynomial. 

We next compute the rational generating function 

5(y Parcto n( . +JjBoo);z) 

from g(V Pareto ; z) using the Intersection Lemma. Let e' > be a rational 
number, which we will determine later. By Theorem 14, we compute a 
solution v e / G y Pareto with 

/(vv)>(l-e')/(v*), 

or, equivalently, 

/(v*) - /(v e ,) < e'/(v*). 

Thus, 

[dg(v, v,)] D - [d Q (v, v*)] D = ||v e , - v-||g - ||v* - v-||g 

= /(V*) - /(VeO 

<e'/(v*) 

= e '((/? 7 /« 2 ) D -K-v||g). 

Since 7 is the smallest integer with (22) and also ||v* — v|| is an integer, 
we have 

7 < ||v* - vl^ </3||v*-v||q. 

Thus, 

2D 

[do(v,v ( .)] D -[<io(v,V)] D < e '' 



l v - V II Q 
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An elementary calculation yields 



d<2(v,v e /) - d Q (v,v*) < — 



2D 



Thus we can choose 



eD 







a 



2D 



-i -1 



1 



(25) 



to get the desired estimate. Since a, (5 and D are fixed constants, we have 
e' = 0(e). Thus the computation of vy S yParcto ^ Theorem 14 runs in 
time polynomial in the input encoding size and -. □ 

Remark 15. It is straightforward to extend this result to also include the l v 
norms for odd integers p, by solving the approximation problem separately 
for all of the 2 k = 0(1) shifted orthants v + 0„ = { v : cr^Uj - vi) > 0}, 
where a G {±l} fc . On each of the orthants, the £ p -norm has a representation 
by a polynomial as required by Theorem 13. 
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